Ten-micron Observations of Nearby Young Stars 
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We present new 10/mi photometry of 21 nearby young stars obtained at the 
r£3 ■ Palomar 5-meter and at the Keck I 10-meter telescopes as part of a program to 

search for dust in the habitable zone of young stars. Thirteen of the stars are 
in the F-K spectral type range ("solar analogs"), 4 have B or A spectral types, 
and 4 have spectral type M. We confirm existing IRAS 12/im and ground-based 
10/iHi photometry for 10 of the stars, and present new insight into this spectral 
regime for the rest. Excess emission at lOyum is not found in any of the young solar 
analogs, except for a possible 2.4-sigma detection in the G5V star HD 88638. The 
G2V star HD 107146, which does not display a 10/xm excess, is identified as a new 
Vega-like candidate, based on our 10/im photospheric detection, combined with 
previously unidentified 60/zm and 100/im IRAS excesses. Among the early-type 
stars, a 10/im excess is detected only in HD 109573A (HR 4796A), confirming 
prior observations; among the M dwarfs, excesses are confirmed in AA Tau, CD - 
40°8434, and Hen 3-600A. A previously suggested N band excess in the M3 dwarf 
CD -33° 7795 is shown to be consistent with photospheric emission. 

We calculate infrared to stellar bolometric luminosity ratios for all stars 
exhibiting mid-infrared excesses, and infer the total mass of orbiting dust in 
the cases of optically thin disks. For a derived median photometric precision 
of ±0.11 mag, we place an upper limit of M dust « 2 x 10~ 5 A<f e on the dust 
mass (assuming a dust temperature of 300 K) around solar analogs not seen 
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in excess at 10/xm. Our calculations for the nearby K1V star HD 17925 show 
that it may have the least massive debris disk known outside our solar system 
(M dust >7x 10~ 6 M e ). 

Our limited data confirm the expected tendency of decreasing fractional dust 
excess fd = Lir/L* with increasing stellar age. However, we argue that esti- 
mates of fd suffer from a degeneracy between the temperature and the amount 
of circumstellar dust M dust , and propose a relation of M dust as a function of age, 
instead. 

Subject headings: circumstellar matter — dust, extinction — planetary systems: 
protoplanetary disks — infrared: stars 

1. INTRODUCTION 

Since the discovery of far-infrared (far-IR) excess emission by the Infrared Astronomy 
Satellite (IRAS) toward the main sequence star Vega (Aumann et al. 1984), nearly 400 
other main-sequence stars have been detected with such excesses (Mannings & Barlow 1998; 
Sylvester & Mannings 2000; Habing et al. 2001; Spangler et al. 2001; Laureijs et al. 2002a; 
Song et al. 2002, and references therein), interpreted as due to orbiting dust. Too old to 
possess remnant primordial dust (which would be cleared by Poynting-Robertson drag on 
the time-scale of several Myr in the absence of gas), these stars owe their far-IR excess 
to emission by "debris disks," formed by the continual collisional fragmentation of larger 
bodies (Backman & Paresce 1993, and references therein). Stellar systems exhibiting this 
phenomenon have been accordingly named "Vega-like." Subsequent higher spatial resolution 
imaging at millimeter wavelengths of the nearest subsample of Vega-like stars has revealed 
intricate disk-like structures, with gaps and concentrations indicative of ongoing disk clearing 
by planets (e.g., Holland et al. 1998; Greaves et al. 1998; Jayawardhana et al. 1998; Koerner 
et al. 1998; Wilner et al. 2002). 

While such observations have provided a wealth of information about the outer reaches 
(> 10 AU) of these systems, the "habitable zone" (defined as the region corresponding to 
the temperature of liquid water) around the majority of them remains at separations un- 
resolvable by direct imaging, and can be probed only by interferometry or by mid-infrared 
(mid-IR) spectroscopy. To first order, however, information about the presence of warm 
(~300 K) material around a star can be extracted from its spectral energy distribution 
(SED), and the presence or absence of a mid-IR excess above the expected photospheric 
level. 
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To date, searches for mid-IR excesses from both the ground and space have been domi- 
nated by upper limits. All Vega-like stars have been identified from IRAS and Infrared Space 
Observatory (ISO) observations, and are therefore limited by their sensitivities (~0.5 Jy at 
12-60/mi for the IRAS Point Source Catalog). In a sample of nearly 150 stars in the 1- 
630 Myr age range, Spangler et al. (2001) find that the vast majority of targets older than 
~20 Myr have either no measurable 12 and 25 /im excesses or no detections from IRAS. Nev- 
ertheless, given a gradually decreasing dust fraction for young stars (Spangler et al. 2001; 
Habing et al. 1999), more sensitive mid-IR observations are expected to uncover less massive 
debris disks around older stars. 

Young stars form a particularly suitable set of targets to search for mid-IR emission 
from circumstellar disks. They are likely to be at a transition stage, where the innermost 
regions of the dust disk have been cleared, whereas, due to an increase of the viscous and 
dynamical timescale with radius (Hollenbach, Yorke, & Johnstone 2000; Nakagawa, Hayashi, 
& Nakazawa 1983), cooler material may still be present farther away. No excess is expected 
at near-infrared (near-IR) wavelengths (1-5/im; < 0.1 AU) from these stars, yet longer- 
wavelength excess due to more distant (0.1-10 AU), cooler material, could be detectable. 

Except for stars within a few tens of parsecs, stellar photospheres are not detected at 
12/zm with the small apertures of IRAS and ISO. However, the recent availability of mid-IR 
instrumentation on 5-10-meter class ground-based telescopes has made them competitive 
with IRAS/ISO within the atmospheric windows at 10/xm and 20/xm, with limiting flux 
densities <20 mJy detectable at 10/im. In particular, for sufficiently nearby and sufficiently 
young stars, sensitivity to stellar photospheres is achieved. This implies sensitivity to a small 
excess above the photosphere, and hence to debris disks ~ 10~ 5 M e , limited by photometric 
precision. In addition, the higher angular resolution attainable with large ground-based 
telescopes allows for unambiguous determination of the origin of the infrared (IR) flux within 
the IRAS/ISO beam. We can thus discriminate between Vega-like excesses associated with 
the star, and apparent excess flux associated, instead, with a background object in the IRAS 
or ISO beam (Aumann & Probst 1991; Lisse et al. 2002), structure in the IR cirrus, or due 
to detection threshold bias (Song et al. 2002). 

In this paper, we present 10/im photometry of 21 nearby stars spanning the ~l-600 Myr 
age range, and thus expected to cover the stages in the transition between a primordial and a 
debris disk. Age estimates are obtained from the literature, and are based on high-resolution 
spectroscopic measurements of lithium abundances, and/or from inferred association with 
known kinematic groups of stars. Our sample was selected to include nearby (< 150 pc) 
stars, covering a wide range of spectral types (B9-M4), for which the photosphere should be 
detectable at 10/zm from the ground. 
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Twelve of our targets have IRAS and/or ISO detections, 8 of which have already been 
inferred to harbor debris disks. For these objects, our 10/mi data test the existence of possible 
excess at shorter wavelengths. 

2. GROUND-BASED TEN-MICRON PHOTOMETRY 

Observations were conducted under marginally photometric conditions at 10.7/xm with 
the Long Wavelength Spectrograph (LWS; Jones & Puetter 1993) on the Keck I 10-meter 
telescope during two runs (2000 February 20-21 and December 9), and under photometric 
conditions at 10.3yUHi with SpectroCam-10 (SC-10; Hayward et al. 1993) on the Hale 5-meter 
telescope, on 2002 January 1 and 3. The 10.7/im LWS filter has a full-width half-maximum 
(FWHM) of ~ 1.6/im 1 , and the 10.3/im SC-10 filter has a bandwidth of ~ 1/im (Hayward, 
Houck, & Miles 1996), and are thus both longward of the telluric ozone feature at 9.6/im, 
and of the peak of the potential photospheric silicate emission at 9.7/xm. 

Table 1 lists the instrument and epoch of observation for each source. 

2.1. LWS Data 

LWS is a Boeing-made 128 2 pix Si:As moderate-flux array, mounted at the Forward 
Cassegrain module on the Keck I telescope, providing diffraction-limited imaging (10" field 
of view, 0.08"/pixel) and spectroscopic (12=100 and 1400) capabilities in the 3.5-25/im wave- 
length range. For the observations presented here we used only the imaging mode of LWS 
at 10.7/im. The seeing was poor and variable: 0.3"-0.6" at 10/xm (the attainable diffraction 
limit is 0.22"). Data were taken in the standard manner employed for mid-IR observations, 
with both chopping and nodding to off-source positions. The resulting images were pro- 
cessed and co-added using the LWSCOADD routine written in IDL by Gregory D. Wirth. 2 
No flat-fielding was performed, as this was found to increase the scatter in the photometry. 
Instead, we attempted to place the stars always in the same place on the array. 

Magnitudes were obtained based on measurements of 10 IRTF iV-band standard stars 
(Tokunaga 1984; Rieke, Lebofsky, & Low 1985). Photometry was performed with the 
IRAF/PHOT task in apertures of radius 0.96" (12 pix), with sky background measured 



1 http://www2. keck, hawaii.edu/inst/lws/filters. html 

2 Available at http://www.astro.caltech.edu/mirror/keck/inst/lws/lwscoadd_pro.html 
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in 2.0"-3.2"-radius (25-40 pix) annuli. A curve-of-growth correction (—0.12 ± 0.01 mag, 
based on standard star measurements) to a 1.92" (24 pix) aperture was then applied to each 
star to bring the magnitudes to within ^2% of the infinite aperture value. An atmospheric 
extinction correction of —0.32 ±0.07 mag/airmass was derived based on standard star obser- 
vations, where we have allowed for ^20% shifts in the night-to-night atmospheric zero-points 
to align the extinction curves from each night. The root-mean-square (r.m.s.) scatter in the 
photometry of the standards is 0.06 mag for 2000 February 20 and 21, and 0.13 mag for 2000 
December 9. 



2.2. SC-10 Data 

SpectroCam-10 is a 128 2 pix Si: As BIBIB (Back Illuminated Blocked Impurity Band) 
detector manufactured by Rockwell. For our program we used it in imaging mode (16" 
field of view) with a 10.3/im filter, although the instrument also allows i?=100 and 2000 
spectroscopy between 5 and 20/xm. The seeing was poor, 0.7" at 10/im, but relatively stable. 
Both chopping and nodding were used in the data acquisition, with the chop and the nod 
throws being of the same magnitude (^20"). The chip was binned into a 64x64 pix array 
with 0.256" pixels, providing nearly Nyquist sampling of the diffraction limit (0.48") of the 
telescope, and then stored into a 3-dimensional image, each slice of which contained a co- 
added chopped frame. Dark frames were taken for all source exposure times, and subtracted 
from all frames. Flat fields were constructed by median-combining the off-source exposures 
for each target. This approach of obtaining a flat field was preferred to making a master 
flat from all exposures during the night, because of an unexplained ~4" shift in the filter 
vignetting pattern across the field of view over the course of the observations. The r.m.s. 
scatter in the flat fields was <1%. 

Final images from SC-10 were obtained by differencing each nodded pair, and then 
averaging the differences. Since source re- acquisition during nodding was not perfect, the 
PSF of the final image was broader (FWHM ps 3 pix = 0.77") than the diffraction limit. 
Attempts with centroiding the target in each differenced pair before co-adding were not 
successful, because only the standard stars had enough signal-to- noise (S/N) to obtain good 
centroids. Object magnitudes were measured on the reduced images in 1.5"-radius (6 pix) 
apertures, using the IRAF/PHOT task. Local sky was estimated as the mode of the pixels 
in a 3.8"-5.1" (15-20 pix) annulus around each object. An aperture correction (—0.147 ± 
0.014 mag, determined from measurements of standard stars) to a 3.1" (12 pix) aperture was 
then applied to bring the estimate of the enclosed flux to within ^2% of its infinite aperture 
value. In the case of the double source RE J0137+18A/B (1.7" separation) magnitudes were 
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obtained through PSF-fitting with DAOPHOT/ALLSTAR (Stetson 1987), using a 3.1"- 
radius PSF fit to the standard star (3 And, observed immediately beforehand. An airmass 
correction of — 0.13±0.03 mag/airmass, as determined from observations of the standard 
stars a Boo and (3 Gem on 2002 January 3, was applied to all measured magnitudes. The 
r.m.s. scatter in the photometry of the standards is 0.015 mag for 2002 January 1, and 
0.007 mag for January 3. 



3. SEDS AND PHOTOSPHERIC MODELS 

3.1. Optical and Infrared SEDs 

SEDs for the observed stars are presented in Figure 1, with our 10//m photometry 
plotted as solid circles. UBVRILMNQ measurements were collected from the literature, 
where available. JHK S photometry was extracted from the 2MASS Point Source Catalog 
for all but the brightest sources (HD 216803, HD 102647, HD 17925, and HD 109573A), the 
data for which were taken from Aumann & Probst (1991) and Jura et al. (1993). For the 
UBVRI and JHKLMNQ sets of filters, conversion from magnitudes to flux densities Fx 
was performed using the zeroth- magnitude flux densities listed in Allen (1973) and Cohen 
et al. (1992), respectively. For the objects with 2MASS photometry, the transformation 
from JHK S magnitudes to flux densities was based on the calibration by Cohen, Wheaton, 
& Megeath (2003), outlined in the 2MASS Explanatory Supplement (Cutri et al. 2003, 
ch.VI.4a). For the N filter we have adopted the same conversion factor both for the 10.3/wn 
(SC-10) and for the 10.7/xm (LWS) filters. Given the difference in central wavelengths of 
the two A^-band filters, the derived zero- magnitude flux densities are expected to be within 
^4% of each another. 

IRAS fluxes and upper limits were obtained from the IRAS Faint Source Catalog 
(FSC), or from the Point Source Catalog (PSC) when the FSC did not contain the object 
(HD 147809, HD 109573A, and the 100//m data point for CD -40°8434). FSC photometry 
was given preference over the PSC because of its higher signal-to-noise, and superior sen- 
sitivity to faint objects. Reduced ISO data in the 12-60/iin range were available for 5 of 
our targets. For HD 17925 and HD 102647 we adopted the color-corrected 25/zm ISO flux 
densities presented in Laureijs et al. (2002a); the 11.5/xm and 20/iin data for HD 216803 are 
from Fajardo-Acosta et al. (1999), and 60/iin measurements for HD 145519 and HD 147809 
were taken from M. Meyer et al., in prep. The ISO Data Centre lists that additional data 
were obtained for 5 of our targets (HD 147809, HD 145519, HD 216803, CD -40°8434, and 
AP 93), none of which were user- reduced, however. For the sources with no available reduced 
IRAS or ISO data we have adopted nominal IRAS all-sky survey detection limits away from 
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confused regions of the sky: 0.5 Jy at 12, 25 and 60/im, and 1.5 Jy at 100//m. 

The IRAS and ISO flux densities were color-corrected assuming blackbody emission, 
following the procedures outlined in the IRAS Explanatory Supplement (Beichman et al. 
1988, ch.VI.C.3) and the ISO Handbook (Laureijs et al. 2002b, vol.IV, Appendix C). The 
adopted color correction factors were chosen depending on the effective temperature T c s of 
each star (based on spectral types from Simbad, unless otherwise noted; Table 2). For IRAS 
they were 1.41-1.45 at 12/im, 1.39-1.41 at 25/im, 1.30-1.32 at 60/im, and 1.09 at 100/im, 
where the range of values corresponds to the range of T eff of the stars in our sample. 

Near- and mid-IR photometry for the observed sources from 2MASS (JHK S ), IRAS 
(12/im), and from our 10/im observations is summarized in Table 2. Twelve-micron color- 
corrected magnitudes ([12]) were calculated assuming an IRAS color-corrected flux of 28.3 Jy 
for a zeroth-magnitude star (Beichman et al. 1988, ch.VI.C.2a). 

3.2. Photospheric and Dust Models 

The stellar SEDs in Figure 1 were fit with photospheric models from Kurucz (1979, 
1992) and with NextGen models from Hauschildt, Allard, & Baron (1999). In fitting model 
photospheres, we have assumed log (7 = 4.5 and solar metallicity in all models, as is ap- 
proximately true for main-sequence dwarfs in the solar neighborhood (Drilling & Landolt 
2000). The temperatures of the models were chosen to correspond to the spectral types of 
the stars (Table 2; de Jager & Nieuwenhuijzen 1987). This produced satisfactory (by eye) 
fits in most cases, except when non-simultaneous photometry of variable stars in different 
filters had to be fitted (e.g., AA Tau, AP 93). For RE J0137+18A/B we found that the 
required photospheric temperature (T eff rs 4700 K for spectral type K3V) was too high to fit 
the optical and near-IR SED. We obtain better fits using a NextGen 4200 K photosphere for 
both stars, which can be explained if both components of the binary had sub-giant gravities 
(logg m 4.0; de Jager & Nieuwenhuijzen 1987). Given their young age (~ 1 Myr), AA Tau, 
and DI Tau are also expected to have sub-giant gravities. However, the effect on the model 
photospheres is small compared to that due to the extinction measured toward both stars 
(A v = 1.0 and 0.76, respectively; Figure lb). 

The agreement between the two sets of models is satisfactory over the 5000-10000 K 
range of temperatures, though at lower temperatures the NextGen models describe the 
stellar photosphere more closely (Hauschildt et al. 1999). For stars with T cff > 5500 K, we 
have chosen to plot Kurucz models, as their lower resolution gives a better visual idea of 
the integrated flux density of the stellar photosphere at short wavelengths. At wavelengths 
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A > 20/xm, for which Kurucz does not predict fluxes, we have linearly extrapolated the model 
photospheres from the last two (10/im and 20/xm) data points. NextGen models (ranging 
from 10 nm to 10 mm) have been used for all cooler (T cff < 5000 K) stars. 

The model for every star was normalized to each of the measured J, H and K fluxes 
(although in Figure 1 we have shown only the J-band normalizations). The mean predicted 
photosphere at 10.3/im (for the SC-10 targets) or at 10.7/im (for the LWS targets) from the 3 
normalizations was then taken as representative of the photospheric level, and the standard 
deviation of the mean - as the error (ajvpred) in the predicted photosphere. The near-IR 
was chosen for normalization because: (1) it is close in wavelength to the target band at 
10/im; (2) uniform photometry is available for the majority of our targets (from 2MASS); 
and (3) most of our targets were not expected a priori to possess near-IR excesses. 

We have assumed no circumstellar extinction except in the cases for which an estimate 
was found in the literature: HD 145519 (A v = 1.0, Krelowski & Strobel 1983), HD 109573A 
{A v = 0.03, Jura 1991), HD 147809 {A v = 1.6, Eggen 1998), AA Tau {A v = 1.0, Beckwith et 
al. 1990), AP 93 (A v = 0.31, O'Dell, Hendry, & Collier Cameron 1994), DI Tau (A v = 0.76, 
Stassun et al. 2001), and CD -40°8434 (A v = 0.74, Gregorio-Hetem & Hetem 2002). An 
extinction law was adopted from Mathis (1990) over the entire SED range (0.25/im-10mm). 

Table 2 lists the derived 10/im excesses AA for all targets. The quoted error in AA 
is comprised in each case of the photometric error in the measurement of A, and the error 
crjvpred in the level of the predicted photosphere. Excluding the stars already known to 
possess such (HD 109573A, CD -40°8434, AA Tau, and Hen 3-600A), and AP 93 (for 
which we only measure an upper limit to the 10/im intensity), the mean A-band excess is 
0.01 ± 0.03 mag. This is indicative of the accuracy with which the photospheric models are 
able to predict the broadband 10/im flux density, and demonstrates that the models do not 
systematically over- or under-estimate the stellar photosphere at A band. 

Excess IR emission, whenever apparent in the SED of a star, was fit by a single blackbody 
curve, as would be expected from a narrow ring of dust at a fixed distance from the star. 
The temperature of the blackbody curve (T^) is identified for each case in Figure 1. The 
fits were performed by eye, and experimentation with a range of values for T bb showed that 
the determined temperatures were generally accurate to ±10%. For CD -40°8434, AA Tau, 
and Hen 3-600A, single temperature blackbody emission could not account for the excess at 
wavelengths >60/iin, and for these stars we have allowed for additional, cooler blackbodies 
(T bbl ,T bb2 , etc). Although a disk model would be more appropriate for fitting the observed 
excesses, our simplistic treatment allows us to estimate color correction factors for the IRAS 
and ISO data (which were then re-applied instead of the ones based on the stellar T e g), and 
to obtain approximate measurements of the dust luminosity fraction in each case. 
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HD 17925 and HD 145519, with significant excesses only at 60yum, are the only stars for 
which we have not applied any color corrections at wavelengths > 60/im. The large uncer- 
tainty in the blackbody temperatures (30-250 K for HD 17925, 30-140 K for HD 145519) 
that fit the far-IR excesses and upper limits for these stars correspond to a range of values for 
the color correction (0.91-1.19 for IRAS at 60/im for HD 17925, 0.93-1.08 for HD 145519). 
Thus the plotted (Figure 1) IRAS and ISO flux densities and upper limits for these stars 
at wavelengths > 60/im correspond to flat spectrum sources (F u oc z^ 1 ), as is the default 
calibration for IRAS and ISO. 



Ten of our sample stars have been previously reported to possess mid- or far-IR dust 
excesses: HD 145519, HD 109573A (HR4796A), HD 147809, HD 102647 ((3 Leo), HD 17925, 
AA Tau, DI Tau, CD -40°8434, Hen 3-600A, and CD -33°7795. Three among these 
(HD 145519, HD102647, and HD 17925) have known 60/im dust emission from ISO and/or 
IRAS (Backman et al. 1997; Laureijs et al. 2002a; Habing et al. 2001), but do not exhibit 
significant excesses at 10/xm, in agreement with previous measurements. A fourth, DI Tau, 
previously suspected of harboring warm circumstellar dust on the basis of 12/im IRAS data 
(Skrutskie et al. 1990), has subsequently been shown not to possess an excess at 10/im (Meyer 
et al. 1997), and that all of the excess IRAS flux can be attributed to the neighboring (15.1") 
DH Tau. Our photometry of DI Tau is consistent with that of Meyer et al. (1997) and with 
the lack of circumstellar dust. 

Of the remaining 6 stars that have been reported previously to possess excesses at lOyum 
or 12/im, we confirm such in 4 of them: HD 109573A (Jura et al. 1998), AA Tau (Strom et 
al. 1989), CD -40°8434 (Gregorio-Hetem et al. 1992), and Hen 3-600 A (Jayawardhana et al. 
1999a). Our A-band data for these 4 objects agree very well with the previously published 
10/iHi photometry, as well as with the color-corrected IRAS 12/im flux densities, and we will 
not discuss their photometry further. 



4. 



RESULTS 



4.1. 



Stars with Confirmed Excesses 



4.2. 



Stars with Unconfirmed Excesses 



There are 2 stars — CD -33°7795, and HD 147809 — for which our non-detection of 
excess at 10/iin is in disagreement with previous data and interpretations. 
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4.2.1. CD -33° 1195 

Jayawardhana et al. (1999b) report a possible 2.6-sigma detection of a 10/zm excess in 
the M3V star CD -33°7795. Their claim is based on a measured flux density of 96±9 mJy, 
and an estimated photospheric level of 70±5 mJy at 10.8/xm, assuming K — N = 0.0. Our 
10.7/im flux density (93±8 mJy) is consistent with theirs, and with the level of photospheric 
emission (100±5 mJy) found from the NextGen model. The predicted K — N color by the 
NextGen model is 0.27, which agrees with the measured K s — N = 0.24 ± 0.04 for this 
star. For comparison, Kenyon & Hartmann (1995) list K — [12] = 0.49 mag for M3V main- 
sequence stars, which is expected to be comparable to their K — N color. There is thus no 
indication of excess emission at 10/im in this object. Jayawardhana et al. acknowledge that 
if K — N ~ 0.3 their measurement would be consistent with being photospheric. 



4.2.2. HD 141809 

Our 10.7/xm measurement of 46±6 mJy for this star, while consistent with being pho- 
tospheric, is a factor of 5 fainter than the IRAS flux density of 243±24 mJy at 12/im. 
HD 147809 is detected at 12/im and 60/im (S/N 5 in both bands), but is listed only 
in the PSC, despite being fainter than the nominal PSC sensitivity at these wavelengths 
(~0.5 Jy), and brighter than the sensitivity limit of the FSC (~0.2 Jy). However, the IRAS 
Faint Source Catalog Rejects (FSCR) lists a 12/im source (IRAS Z16219-2514) within 5" of 
HD 147809 (associated with IRAS 16219-2514): a separation equal to a fifth of the angular 
resolution of IRAS at 12/im. The reason for rejection (and inclusion in the FSCR) is that 
IRAS Z16219-2514 is detected only in 1 band (12/im) in the FSC, with signal-to-noise of 3-6 
(Moshir et al. 1992, actual S/N = 5.9). The FSCR also lists one nearby 100/im-only cirrus 
extraction ("cirrus" flag 1), which may have possibly affected the PSC 60/xm measurement. 
The possibility that the 60/xm and the 12/xm PSC detections are associated with IR cirrus 
is explored below. 

HD 147809 is a member of the Upper Scorpius OB complex in the p Oph region, which 
is associated with highly- variable background emission in the mid- and far-IR (Ryter, Puget, 
& Perault 1987), as seen in the IRAS maps of this area. Backman et al. (1997) observe the 
star with ISO at 60/im and 90/rni with a detection only at 60/im: 0.91±0.07 Jy (S/N— 13) - 
4-sigma fainter than the 1.64 ±0.18 Jy flux listed in the IRAS PSC at this wavelength. This 
is possible if the peak of the 60/im source was in the central (45x45") pixel of the 3x3 
PHT-C100 ISO array, and if the emission extended over the 8 border pixels, which were 
used to estimate local background. IRAS photometry, on the other hand, does not involve 
local background subtraction (each survey scan is initiated and terminated with a flash of 
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the internal reference source to monitor responsivity of the system; Beichman et al. 1988, 
ch.VI.B.l.), and thus measures the total flux in the IRAS beam (60" at 60pm; Beichman et 
al. 1988, ch.II.C.3.). The lower ISO flux density is therefore indicative of the extended nature 
(over several C100 pixels: >50") of the 60pm emission. At the distance of the Upper Sco OB 
association (~160 pc; de Geus, de Zeeuw, & Lub 1989; Jones 1970), the physical size of the 
structure would be ~10000 AU. While it is possible that HD 147809 possesses a dust disk of 
such size, the corresponding dust temperature (10-30 K) is inconsistent with the 90pm ISO 
upper limit (read off from the default ISO post-stamp measurement representation 3 ) 

If the 12pm excess in HD 147809 were associated with orbiting circumstellar mate- 
rial around the star, the ratio of the PSC 12pm and 25pm (upper limit) flux densities 
S*,, (1 2pm) /S^ (25pm) > 0.95 ± 0.09 would require a blackbody temperature T bb > 154 K 
(3-sigma lower limit), whereas the measured S u (l0.7fim)/S u (l2fxm) = 0.19 ± 0.03 requires 
T b b = 65 ± 5 K. Therefore, the 12pm emission arises either outside our l"-radius LWS aper- 
ture, or is extended over at least several arc seconds, s.t. its intensity in the LWS aperture 
is comparable to that in the background annulus (2.0-3.2"). An S u (12fim)/S u (25nm) ratio 
greater than unity may be indicative of non-thermal emission from aromatic infrared lines 
in the 12/im band, as previously suggested in association with inter-stellar cirrus in pOph 
(Meyer, Wilking, & Zinnecker 1993; Ryter et al. 1987). 

Summarizing the evidence presented above, we conclude that the 12/im and 60/im ex- 
cesses in the SED of HD 147809 are not associated with the star, but are likely due to 
inter-stellar cirrus. 



4.3. New Excess Detections 

We find previously unknown IR excesses in 3 of our target stars. HD 88638 is the only 
one for which we present new evidence for possible excess emission at 10pm. HD 107146 is an 
overlooked Vega-like star from IRAS, and HD 70516 has a 29" distant companion probably 
responsible for the 12pm IRAS excess. 

4.3.1. HD 88638 

The flux measured from this star at 10.3/im is 2.4-sigma (AiV = 0.19 ± 0.08 mag) 
above the 5500 K, solar metallicity and gravity (log (7 = 4.5) Kurucz photosphere. Given 



3 Available at the ESA ISO website, http://www.iso.vilspa.esa.es/ida/index.html 
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the Hauschildt et al. (1999) claim that the NextGen models better approximate the stellar 
photosphere at effective temperatures <5500 K, and the spectroscopically determined T eS = 
5360 K for this star (Strassmeier et al. 2000), we also fitted a 5400 K solar gravity and 
metallicity NextGen model to the SED. The observed excess is 2.3-sigma over the NextGen 
model: AN — 0.18 ± 0.08 mag; i.e. the excess is nearly model-independent. 

There is no measurement of the gravity or metallicity of HD 88638 in the literature. 
Given its Hipparcos distance (38 ±4 pc) and magnitude (V = 8.02), and a bolometric correc- 
tion of -0.165 (Flower 1996), the star should have a luminosity of = 0.83 ± O.15L , and a 
radius of i?* = 1.07 ± 0.10i? Q . Assuming ~ 1M for a G5 star of this luminosity, its gravity 
is log g = 4.4, i.e., nearly solar. Alternatively, if sub-solar metallicity were causing the higher 
10/im flux, it would have to be < —2.0 dex for the Kurucz model photosphere to fall within 
lei of our N band measurement. Such low metallicity would imply that HD 88638 is a halo 
(population II) star, whereas its space velocity (U,V,W) = —(62.6,8.5,-1.3; Strassmeier 
et al. 2000) places its orbit in the plane of the galactic disk (population I). Gravity and 
metallicity effects are thus unlikely to be the cause for the observed excess. 

Interstellar extinction of order A v ps 0.1 could contribute to the abnormal mid-IR 
emission. However, HD 88638 is well out of the galactic plane (b = +53°46'), and at 38 pc 
from the Sun extinction along the line of sight is expected to be negligible. Moreover, 
reddened photospheric models normalized to the near-IR tend to underestimate the optical 
flux by ~0.2 mag. Hence, extinction is not the probable cause. 

The blackbody temperature range that accommodates the excess, the IRAS upper limits, 
and the near-IR data is 130-1500 K. This does not exclude a substellar companion of spectral 
type L6 or later (T e s < 1500 K). For an L6 dwarf Mk > 11.4, and assuming K—N ps K—L' > 
1.3 (Leggett et al. 2002), at the distance of the primary (distance modulus 1.4 mag) the 
potential companion would have N > 11.5. The 10/im excess that this would produce in 
HD 88638 (N = 6.04, Table 2) is AN < 0.01: far too faint to account for what is observed. 

Therefore, if real, the 2.4-sigma N band excess is most likely due to orbiting debris. 
However, from a statistical point of view, a 2.4-sigma positive deviation should occur in 
0.82% of all observations, if the errors are Gaussian. Given 21 observed sources, a random 
2.4-sigma excess will occur in 16% of such sets. There is a hence non-negligible probability 
that the observed 10.3/zm excess is random. 

Due to the marginality of our detection, and to the large range of possible blackbody 
temperatures, we have not re-applied the corresponding color-correction factors to the IRAS 
upper limits. These would range from 0.90-1.36 at 25/im and from 1.08-1.29 at 60/zm (the 
ones applied based on a 5500 K stellar SED are 1.40 and 1.32, respectively), and will have 
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the effect of allowing a ~10 K cooler blackbody fit, i.e., 120 K. 

4.3.2. HD 107146 

Although we do not detect above-photospheric emission at 10/zm from this source, we 
choose to discuss it here because it has not been previously identified as a Vega-like star, 
despite 60fim and 100/xm excesses measured by IRAS. Because its 12/im IRAS flux is con- 
sistent with being photospheric, and because of a non-detection at 25yum, it would have 
been omitted in most previous searches concentrating on the 12/xm and 25/zm IRAS bands. 
However, this object has also been missed in the Stencel & Backman (1991) survey for IR 
excess emission from IRAS 12/j.m point sources at high galactic latitude (|6| > 25°). This 
survey selected stars with 12/xm IRAS detections from the PSC that were also detected in 
at least one more IRAS band, and that had positional association with an SAO star. It was 
apparently tailored to discover excess IR emission in sources like HD 107146 (SAO 100038, 
6=77.04), detected in the 12/im, and in the 60/im bands in the PSC. However, the low signal- 
to-noise ratios of these detections, S/N=3.9 and 6.9, respectively, provide a clue as to why 
the star may have been overlooked. No ISO data exist for this source. 

The detection of an excess in this star is not false, as has been suggested for other re- 
ported candidate Vega-like stars with \ow-S/N excess in a single IRAS band. The misiden- 
tification in those cases is explained as a likely result of detection threshold bias (Song et al. 
2002). In the case of HD 107146 however, the FSC (which contains higher signal-to-noise 
measurements for objects fainter than the PSC completeness limits) lists detections of the 
star in 3, rather than 2 bands. The higher significance of the detection at 60/im (S/N =20), 
and a 100/im (S/N =6) detection make it highly improbable that the excess is due to an 
upwards noise fluctuation. Moreover, our ground-based measurement is consistent with the 
listed IRAS 12//m flux density, and with being photospheric, thus excluding the possibility 
that the IRAS flux comes from an unassociated background object in the IRAS aperture. 

The IR excess in this star can be fit by a 60±10 K blackbody. See Williams et al. (2003) 
for a further discussion of this source, including recent sub-millimeter observations. 

4.3.3. HD 70516 

The 12/zm IRAS detection of this star is listed with a 10-sigma significance in the FSC, 
and is 4.1-sigma above the NextGen photosphere (A[12] = 0.45 ±0.11), whereas our 10.3/im 
measurement is photospheric (AN = 0.07 ±0.09). HD 70516 (HIP 41184; GO) is listed as a 
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proper motion binary in the Hipparcos catalog, with a Ay = 1.5 KO companion (HIP 41181) 
29" away. The separation between the two components matches the resolution of IRAS at 
12/zm (~0.5'; Beichman et al. 1988, ch.I.A.l.), and it is therefore likely that light from both 
components contribute to the measured 12/zm flux density. Given V magnitudes of 7.7 and 
9.2 for the two stars, and expected V — [12] colors of 1.5 and 2.0, respectively (Waters, Cote, 
& Aumann 1987), the combined 12/zm magnitude is 5.85, corresponding to A[12] = 0.35. 
Although a significant fraction of the companion's light will fall outside the IRAS beam at 
12/zm, its vicinity and brightness are marginally consistent with the observed 12/im excess, 
making a debris disk interpretation less likely. 

The binary is not detected at longer wavelengths with IRAS, because its expected 
combined flux density is below the quoted sensitivity limits. 

4.4. Null Excess Detections 

Six of the remaining 8 stars do not have 12/zm measurements from IRAS: RE J0137+18A/B, 
HD 60737, HD 70573, CD -38°6968, and AP 93. The first 5 all have lO^m fluxes consistent 
with being photospheric. AP 93 (observed with the Palomar 5-meter telescope) was not 
detected at 10/zm. The 3-sigma upper limit of 18.8 mJy for this source is 11.0 mJy above 
the model atmosphere, and we can thus exclude any excess emission at this level. 

HD 77407 and HD 216803 have 12/xm measurements from IRAS. Our 10/zm photometry 
for them is consistent with the IRAS data, and the models show that the 10-12/iin emission 
is at the expected photospheric level. 

5. DUST PROPERTIES OF STARS WITH EXCESSES 

5.1. Fractional Dust Luminosity f d 

Table 3 lists the observed dust luminosity fraction f d = L m / L* for all stars with inferred 
disks (marked with "Yes" in the last column). These were obtained based on the fitted 
blackbody curves to the IR excesses. HD 88638 and HD 17925 allow a range of temperatures 
in the blackbody fits, and for them we have listed a range of values for f d , with hotter 
blackbodies generally producing higher fractional luminosities. The measured values for f d 
correspond well to the ones found in the literature (i^(Lit); see the last column of Table 3 for 
references). We are therefore confident that our blackbody fits to the excess IR emission have 
not introduced significant systematic errors in the estimate of fd, with respect to estimates 
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based on more elaborate dust and/or debris disk models (e.g., as in Chiang et al. 2001; Li 
k Lunine 2003). 

For the sources without excess emission at 10/iin we estimate f d < 10~ 3 , approximately 
corresponding to the 3-sigma upper limit on the dust excess that would be produced by 
a hypothesized 300 K debris ring around a G5 star, assuming our median uncertainty of 
0.11 mag in the measured excess AN. This is appropriate, since most of the stars that 
we do not detect in excess are G and K dwarfs. DI Tau is the only M dwarf without a 
detected excess, and given its later spectral type, the upper limit on fd is somewhat higher: 
/d < 3 x 10~ 3 , whereas for HD 147809, the only early-type star (A1V0 not seen in excess, 
fd < 5 x 10~ 4 . For HD 216803, for which the SED is further constrained by the 12/im and 
25/im IRAS photospheric detections, this upper limit is also f d < 5 x 10~ 4 , whereas for 
AP 93, for which we only have an iV-band upper limit, fd < 7 x 10~ 3 . 



5.2. Dust Disk Parameters 

In estimating the mass of the observed circumstellar dust disks, we have followed the 
formalism and the steps outlined in Chen & Jura (2001): 



3L*Q 



a > a min = — (1) 

n(a)da = riQa~ p da (2) 
1 Qty 

M dust > —f dPs D 2 (a) (3) 

t PR = ^^! (4) 

M PB > ^t agc . (5) 

tpR 

Equation (1) gives the minimum particle size a m j n that is not blown out by radiation 
pressure (Artymowicz 1988); and M* are the stellar luminosity and mass, Q pv is the 
radiation pressure coupling coefficient, and p s = 2.5 g cm -3 is the mean density of a dust 
grain (as assumed in Chen & Jura 2001). Equation (2) represents the assumption that the 
particle size distribution n(a)da results from equilibrium between production and destruction 
of objects through collisions (Greenberg & Nolan 1989), where we have assumed p = 3.5 
(Binzel, Hanner, & Steel 2000). If we weight by the number of particles, the equation can 
be integrated to obtain an average particle size (a) = 5/3a min . The inferred dust mass 
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Md us t around the star is found from Equation (3), where we have assumed that the dust is 
distributed in a thin shell at a distance D from the star (Jura et al. 1995). If the grains 
are larger than (a), or if colder dust is present (< 30 K; undetectable by IRAS/ISO), our 
estimate of Md ust is a lower bound. The Poynting-Robertson lifetime tpR of circumstellar 
grains is given by Equation (4) (Burns, Lamy, & Soter 1979). Since this timescale is generally 
shorter than the ages (t age ) of the stars in our sample (Tables 2 and 3), collisions between 
larger bodies are required to replenish the grains destroyed by PR drag. The mass M PB of 
the parent bodies is given by Equation (5). 

Stellar parameters (radius i?*, mass M*, and bolometric luminosity L*) for stars of 
spectral type Kl and earlier were estimated using the Allende Prieto & Lambert (1999) 
database of fundamental parameters of nearby stars, given T e g (based on spectral type) and 
assuming log g = 4.5. For the cooler (T eff < 5000 K) stars these values were taken from the 
literature, or, in the case of Hen 3-600A, and i?* were estimated using the star's distance, 
its visual magnitude (assuming Ay = 0.0), and bolometric correction from Flower (1996). 
Table 5 lists the parameters adopted for each star. 

The radiation pressure coupling coefficient Q pr , which enters in the grain size calculation 
(Equation 1), was assumed ~1, since it deviates from unity only for grains smaller than \ p /2n 
(Burns et al. 1979), where \ p is the wavelength of peak emission from the host star. For the 
coolest K dwarf in our sample (RE J0137+18A/B, assumed T e g = 4200 K) X p /2tt = O.lyum, 
which is comparable to the blow-out size for this star. Our assumption that Q pr ~ 1 is 
therefore valid for K3 and earlier-type stars. 

In deriving the distance D of the dust shell from the star we have expanded the discussion 
in Chen & Jura (2001) by adding a differentiation between grains with blackbody (BB; 
perfect emitters and absorbers), greybody (GB; imperfect emitters, perfect absorbers), and 
inter-stellar medium (ISM; imperfect emitters and absorbers) properties. The three cases 
have been applied as detailed in Backman & Paresce (1993), and we have inverted their 
formulae (3), (5), (6) to solve for L>bb,-Dgb, and -Dism (in units of AU), respectively: 




2 



1/2 



(6) 





In the relations above T dust is the temperature of the dust grains (corresponding to the 
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temperatures T hh of the excesses in Figure 1), and A ~ (a) is the wavelength shortward of 
which the radiative efficiency of a grain of size (a) is roughly constant: e ~ (1 — albedo) ~ 1 
(Backman et al. 1997). Note however, that in fitting the excess in each star, we have only 
used the blackbody law, without modifications at longer wavelengths where particles of size 
a < A emit inefficiently. The effect of this is evident in the SEDs of HD 109573A and AA Tau, 
where at millimeter wavelengths the fitted blackbody overestimates the measured flux. 

The inferred dust properties of the detected circumstellar disks are presented in Table 4. 

It is important to stress here that, given insufficient constraints from the SED, there 
exists a degeneracy between the temperature (Tdust) an d the amount of material (Md ust ) 
around a star required to produce a certain fractional excess fd- This is best seen in the 
case of HD 17925 in Figure 2. The amount of 250 K debris (M dust = 7 x l(T 6 M e , assuming 
GB particles) required to produce the upper limit on the fractional dust excess in this star 
(fd — 2 x 1CT 4 ) is less than what would be needed (Md us t = 4 x 10~ 5 M e ) to produce a 
smaller excess (fd — 6 x 1CT 5 ) at a cooler temperature (70 K). Yet, the two blackbody 
temperatures fit the excess in HD 17925 comparatively well (x 2 = 2.3 and 2.5 per degree of 
freedom, respectively). The best fit (Xp.d.f. = 1-0) * s obtained for a 140 K dust producing a 
fractional excess f d = 7 x 10~ 5 , implying M dus t = 2 x 10~ 6 M e . 

The same trend is observed in the fit to the SED of HD 88638, as discussed in Sec- 
tion 6.1.3 below. We will revisit this observation in Section 6.3 when considering the depen- 
dence of f d and M dust on stellar age. 

6. DISCUSSION 
6.1. Minimum Dust Masses of Stars with Excesses 

6.1.1. Known Optically Thick Disks: CD -4OP84S4, AA Tau, and Hen 3-600A 

The three coolest stars with disks, CD -40°8434, AA Tau, and Hen 3-600A are young 
(< 10 Myr), have high disk luminosity fractions (fd ~ 0.3), and likely possess optically thick 
disks. All of them are, in fact, known T Tauri stars. Our estimated disk radii D therefore 
are not physically meaningful. Detailed modeling of T Tauri disks is beyond the scope of 
this paper, and we have not therefore attempted to estimate grain sizes and disk masses for 
these three stars. For in-depth discussion of each case, see Gregorio-Hetem & Hetem (2002, 
CD -40°8434), Chiang et al. (2001, AA Tau; see also Bouvier et al. 1999), and Muzerolle 
et al. (2000, Hen 3-600A). An estimate of the total disk mass (dust+gas) exists only for 
AA Tau (Beckwith et al. 1990; Chiang et al. 2001), and we have listed that in Table 4. 
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For the remainder of the stars we have assumed that the disks are optically thin, as 
expected given inferred dust luminosity fractions fd < 0.01, higher host star temperatures 
(T eff > 5000 K), and older ages (> 10 Myr). 

6.1.2. Known Optically Thin Disks: ED 145519, ED 109573A, ED 102647, and 

ED 17925 

Four stars have already been inferred to harbor debris disks, three of which (HD 109573A, 
HD 102647, and HD 17925) have been modeled in the previous literature. For these we find 
that our estimates of the disk parameters are in agreement with the existing ones. 

HD 145519. HD 145519 (along with HD 147809) has been included in the analysis of 
Backman et al. (1997) of ISO excesses in open cluster A stars. The star possesses only a 
60/im (3.0-sigma) excess and, unlike around HD 147809, there are no IRAS reject sources or 
cirrus extractions within a 5'-radius. The ISO excess is therefore probably real. The star is 
not discussed in detail in Backman et al. (1997), and no other measurements of the excess 
are found in the existing literature. 

HD 109573A (HR 4796A). Based solely on the aperture size that we have used for N 
band photometry (l"-radius for LWS), we can conclude that the detected 10/zm excess in 
HD 109573A arises within 70 AU of the host star. This is consistent with the extent of 
the disk predicted by the blackbody and greybody particle assumptions, and excludes the 
possibility for the existence of smaller dust particles with properties similar to those of ISM 
grains. Our estimate of the size of the debris disk around HD 109573A agrees with previous 
mid-IR observations (Jayawardhana et al. 1998; Koerner et al. 1998; Telesco et al. 2000, disk 
radii 40 - 70 AU), and with the NICMOS scattered light images of Schneider et al. (1999, 
disk intensity peaking at 70 AU). The dust mass derived here (Md ust > 0.83M ffi = 5.0 x 10 27 g 
for GB particles) is consistent with the one derived in Jura et al. (1995, 3 x 10 26 g < M<j us t < 
8 x 10 27 g), and with the estimates of Li & Lunine (2003, M dust « 4.0 - 7.5 x 10 27 g) found 
from modeling of the near-IR scattered light and of the mid-IR to submillimeter emission 
from the HD 109573A disk. 

HD 102647 ((3 Leo). Our disk parameter estimates for HD 102647 (T dust = 100 ± 10 K, 
Miust ^ 1-7 x 10 _3 M e ) agree with those found in the literature. Laureijs et al. (2002a) 
derive a dust temperature of 83 ± 5 K for HD 102647 from the 25/!m/60/im ISO flux ratio, 
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while Jayawardhana et al. (2001) find T dus t = 120 K from blackbody fits to the photosphere 
and to the IR excess in this star. The amount of dust around HD 102647 inferred from the 
observations of Laureijs et al. is 1.6 x 10~ 3 M ffi : a number not given in their paper, but 
calculated here based on their assumptions (greybody particles) and their formula (A. 5). 
From submillimeter upper limits, on the other hand, Holmes et al. (2003) find Md us t < 
0.59M Mo on = 7.3 x HT 3 M e (M Mo on = M e /81.301; Heiken, Vaniman, & French 1991). Both 
results are consistent with the one presented here. 

HD 17925. Habing et al. (2001) find evidence for a significant ISO excess in HD 17925 only 
at 60/im (with an upper limit on the flux at 90/im), and, assuming a disk distance of 50 AU 
(T eff a: 30 K) from the star (similar to that in the Vega system), find M dust = 1.08 x 10~ 3 M ffi : 
a value which agrees with our range of estimates. However, by considering the full range of 
blackbody temperatures consistent with the IR data, we obtain a range of possible masses 
of the dust disk. Indeed, our measured 1.3-sigma N band excess in HD 17925, combined 
with a 1.3-sigma excess at 12fim (A[12] = 0.09 ± 0.07, from IRAS), and a 2.7-sigma excess 
at 25yum (A[25] = 0.27 ± 0.10 from the averaged ISO and IRAS flux densities, as adopted 
by Laureijs et al. 2002a), point to the possible presence of warmer (140 K) material orbiting 
the star. Recent mid-IR spectroscopic observations by Gaidos & Koresko (2003), also point 
to the marginal presence of a silicate emission peak at 7% of the photosphere at 10/im. 

The IR SED of this star can be fit by a blackbody as warm as ~250 K (Figure 2), 
producing an excess of fd ~ 3 x 10~ 4 . The debris would orbit at a distance of ~ 7 AU if the 
particles are greybodies (~0.8 AU if blackbodies), and would have a mass M dust > 7x 10~ 6 M e 
(lower end of the range listed in Table 4), which could make this the least massive debris 
disk detected outside our Solar System. The required mass of parent bodies would be 
Mp B > 0.015M e , or >3 times the mass of the main asteroid belt in the Solar System (Binzel 
et al. 2000). 

Unfortunately, at the level of the 10-25/xm excesses seen in HD 17925 (l-2cr), we are 
strongly affected by the accuracy of our assumptions for the model photosphere (mainly the 
stellar metallicity), as well as by the precision of the JHK photometry for the object, used 
for the normalization of the model photosphere. Whereas the level of the 10-25/xm excess is 
2.0-2.5o" if the NextGen photosphere is normalized to the J or H apparent magnitude, its 
significance decreases to 0.3-2. Oct if the normalization is done at K. 

A modest amount of extinction (Ay 0.2 mag), a sub-solar metallicity ([Fe/H] < —1.5), 
a sub-dwarf gravity, or a cooler photosphere (by ~ 200 K) could all account for the higher 
mid-IR fluxes in HD 17925. However, given /,j<2x 10~ 4 (Table 3) and the vicinity of the 
star to the Sun (10 ±0.1 pc) we expect Ay < 0.01. The metal abundance and gravity of the 
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star is very nearly solar (Haywood 2001; Cayrel de Strobel & Cayrel 1989; Perrin 1983), and 
its temperature is estimated from high-resolution spectroscopy at 5091 K (Cayrel de Strobel 
& Cayrel 1989; Perrin 1983). None of these factors are thus expected to change the expected 
stellar SED of HD 17925 between 10-25/xm. 

Higher-precision (< 0.10 mag) mid-IR photometry could help constrain the temperature 
and extent of the debris disk around HD 17925, but would be difficult to obtain. However, 
submillimeter imaging can further constrain the SED of the star, and could potentially even 
resolve the emission around it, given its vicinity (10 pc). 

6.1.3. Newly-reported Debris Disks: HD 107146 and HD 88638 

HD 107146. From our analysis of the IRAS FSC data on HD 107146, we report a new 
excess around this star (although, see Williams et al. (2003) for an independent discussion 
of a newly-observed sub-millimeter excess in this source) at wavelengths > 60/im, likely due 
to an orbiting debris disk. Despite a relatively high dust luminosity fraction, fd = 1.5 x 10 -3 
(the IRAS excesses from known Vega-like stars range from 10~ 3 — 10~ 5 ), it has remained 
overlooked due to its apparent faintness. The inferred minimum dust mass Md ust (for GB 
dust particles) is comparable to that in HD 109573A, and should be detectable in the sub- 
millimeter, or in visible scattered light observations. These can then be used to further 
constrain the spatial extent, morphology, and mass of the debris disk. HD 107146 is one 
of very few main-sequence G dwarfs known to harbor debris disks (see Mannings & Barlow 
1998; Sylvester & Mannings 2000; Decin et al. 2000, for other candidates in addition to 
e Eri). At its intermediate age (50-300 Myr; Table 2) it helps bridge the gap between the 
epoch of formation of the Solar System, and its current evolved state of low dust content. 

HD 88638. The evidence for a debris disk around HD 88638 is only tentative, expressed 
as a single 2.4-sigma excess at 10.3/im. The photosphere of the star is below the sensitivity 
limits of the IRAS All-Sky Survey, and it has not been targeted with ISO. Longer- wavelength, 
sensitive IR and sub-millimeter observations (e.g., with SIRTF) are needed to confirm the 
existence of a disk around this star. At the age of HD 88638 (~600 Myr, Table 2), a 
debris disk of the inferred mass and distance from the star (Table 4) would bear the closest 
resemblance to the main asteroid belt in our Solar System among all other known stars with 
debris disks. 

In fitting the excess in the SED of HD 88638, we observe the same trend as discussed in 
the case of HD 17925 (Section 5.2). Because the temperature of the blackbody required to 
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fit the excess is not well-constrained, the disk luminosity fraction f d and the required dust 
mass M dust are not correlated, and are temperature-dependent. 

6.2. Dust Mass Limits for Stars Not Detected in Excess 

With the exception of HD 147809, all the stars in our sample that lack an excesses 
are of spectral type GO and later. The derived upper limits for f d (assuming 300 K debris 
disks) indicate M dust < 10~ 5 M e , assuming that the dust particles behave like greybodies. 
For comparison, the minimum dust mass detectable around a GO dwarf in the ISO 25/im 
survey of Laureijs et al. (2002a) is estimated at 2 x 10 _6 M e , assuming a 120 K disk. Our 
survey appears overall slightly less sensitive than theirs, possibly because of our slightly larger 
median errorbars: ±0.11 mag in AN for our ground-based observations vs. ±0.10 mag in 
A [25] for the ISO observations of Laureijs et al. (2002a). However, if we apply the method 
of dust mass estimation of Laureijs et al., we obtain that our N band observations should 
be sensitive to disk masses M dust > 1 x 10~ 6 M e around GO stars, i.e. our sensitivity is 
comparable to theirs in terms of the minimum amount of detectable dust mass. Regardless, 
from the standpoint of searching for debris disks of different temperatures (and thus, radii), 
the two surveys are complementary. 

6.3. Fractional Excess and Debris Disk Mass vs. Stellar Age 

The results of two large studies targeting disk evolution around main-sequence stars 
with ISO and IRAS data have shown that the amount of circumstellar dust, and the number 
of stars possessing detectable IR dust excesses, are both decreasing functions of stellar age 
(Spangler et al. 2001; Habing et al. 1999, see also Habing et al. 2001). Given our very 
limited sample, and a bias in the source selection toward stars with known dust disks we 
cannot test the conclusions of Habing et al. and Spangler et al. We do, however, reproduce 
the tendency of decreasing dust fraction fd with age among stars with detected IR excesses 
(Figure 3; upper panel), with a power law index of —2.0 ± 0.5 (Xp.d.f. = l-6)> compared with 
-1.76 derived in Spangler et al. (2001). AA Tau and HD 88638 have been excluded from the 
fit, the dust disk of AA Tau being optically thick (and hence not relevant to the discussion of 
debris disk evolution), and the excess in HD 88638 being uncertain. Because we are biased 
towards detecting only the most luminous debris disks, our estimates of fd delineate only 
an upper envelope to the actual amount of IR excess present in any object in the 1 Myr - 
1 Gyr age range. 
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As discussed in Section 5.2 above, the amount of IR excess f d is not in a one-to-one cor- 
respondence with the mass of emitting circumstellar dust M dust , but is also a function of the 
dust temperature, or equivalently, of the distance D of the dust from the star (Equation 3). 
Hotter dust produces larger values of f d by virtue of a larger area under the Planckian 
curve, for a fixed dust mass and mean particle size. It would therefore be a coincidence if 
the /d-age and M dust -age relations followed the same functional form. This would require 
that circumstellar dust is maintained at the same temperature around each star, whereas 
we have examples from our own Solar System (zodiacal dust vs. asteroid belt vs. Kuiper 
belt), and from other Vega-like systems (e.g., ( Lep vs. a Lyr vs. e Eri) that dust temper- 
atures range from ~ 30 — 260 K. One could therefore expect a tighter correlation in the 
log(M dust ) — log(age) vs. the log(fd) — log(age) relation, given that the former is a more 
physical relationship. 

A plot of M dust as a function of stellar age is presented in the bottom panel of Figure 3. 
The fitted slope is log(M dust ) oc (-3.5 ± 0.9) x log(age) (excluding AA Tau and HD 88638), 
with a goodness of fit Xpdf = l-7> similar to the one found for the /d-age relation above. 
There is thus no evidence from our limited amount of data that M dust is more tightly cor- 
related with stellar age than f d . The M dust -age relation is however steeper than the one 
derived by Spangler et al. (2001), and marginally (1-7 a) inconsistent with that expected 
for collisionally replenished secondary dusts disks (power-law index near -2; Zuckerman & 
Becklin 1993). The inclusion of the data point due to HD 88638 only strengthens this trend 
(Xp.d.f. = 1-3 vs - Xp.d.f. = I-? f° r the M dust -age vs. /d-age case), with the power-law indices re- 
maining approximately unchanged (—3.6 ±0.7 and —1.7 ±0.4, respectively), though slightly 
more discrepant with each other (at the 2.4a level, vs. 1.5a when HD 88638 is excluded). 

Given that M dust is a direct estimate of the dust disk mass, and the marginal inconsis- 
tency between the two power law indices, it appears that fd may be a less adequate tracer 
of the evolution of dust disks than our estimate of M dust . The 



M(emitting dust) oc f d oc — — oc t ag 2 (9) 

relation is inferred by Spangler et al. (2001) based on the purely geometric assumption that 
the collisional rate of large particles dN/dt goes as A^ 2 , where N is the number density of 
large planetesimals in the system. When the gravitational potentials Gm/a of the plan- 
etesimals become comparable to the square of their random velocities v 2 ms , the collisional 
cross-sections of the planetesimals increases beyond their simple geometric cross-sections be- 
cause of gravitational focusing. In the two-body approximation, the collision rate is given 
by (e.g., Safronov 1969): 
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2 / 2Gm\ 

a g v rms = 7ia IH — I v rms , (10) 

V fl "mis / 

where <r g is the gravitational cross-section. Collisions between planetesimals become more 
likely with increasing particle size, and their rate goes as a 4 , instead of a 2 as in the geometrical 
regime. The resulting runaway and/or oligarchic growth (Greenberg et al. 1978; Kokubo & 
Ida 1995) should effectively increase the depletion rate of dust particles faster than dN/dt oc 
Cg e ; i-e., 



dN 

M(emitting dust) oc — — oc Cf,, (11) 

dt 



where x > 2. Hence, the steeper power law index that we infer (x ~ —3.5) presents tentative 
evidence for planetesimal growth in debris disks. 

An important shortcoming of the above discussion is that we do not consider the stars 
for which we do not find evidence for orbiting debris disks. As noted, due to an inherent bias 
toward stars with known debris disks in our sample, we cannot perform an analysis similar 
to the ones in Habing et al. (2001) and Spangler et al. (2001). It is significant however, 
that a number of our sources do not show any IR excess, when they would be expected to 
possess such, given universal application of the inferred log(M dust ) — log(age) relation. The 
apparent scatter in disk mass vs. age may depend on the history of each system, and in 
particular, on whether favorable conditions for disk retainment had been present throughout 
the evolution of the star. For instance, the youngest star in our sample, DI Tau, shows 
no excess emission at lOfim, a rather atypical state for an M star < 1 Myr old. DI Tau is 
however a binary system with a ~ 20 AU projected separation, and Meyer et al. (1997) argue 
that the formation of the companion has led to the rapid evolution of the circumstellar disk 
that may have surrounded DI Tau. On the other hand, HD 145519 (Mathieu, Walter, & 
Myers 1989) and HD 109573A also reside in binary systems, though exhibit IR excesses. The 
separation between the components, or conditions other than binarity may be important for 
the disk dissipation time scale. The spectral type of the star may play a role (Decin et al. 
2000), with hotter stars eroding their disks faster, as a result of higher radiation pressures, 
and shorter Poynting-Robertson time scales. In addition, the presence or absence of orbiting 
planets may also affect disk-clearing (Goldreich & Ward 1973). 

A larger, unbiased sample, including a range of stellar spectral types, as well as single 
and multiple systems would be needed to address the issues raised here. More accurate dust 
mass estimates from disk models, employing realistic grains, would be needed to verify the 
validity of our suggestion that M dust (as calculated from Equation 3) is a steeper function of 
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stellar age than fa, and potentially a more adequate tracer of the evolution of debris disks. 

7. CONCLUSION 

Using ground-based 10/zm observations we have demonstrated sensitivity to the stellar 
photosphere of nearby solar-type stars, and hence to small excess (f d = L IR /L* « 1CT 3 
for a 300 K debris disk) above it. Our mid-IR data confirm previous space- and ground- 
based observations for a subset of our targets, and give new insight into this regime for the 
remainder, by complementing the existing data. By employing simple fitting of the observed 
IR excesses with blackbody curves, we extract reliable estimates of the ratio fd of the dust 
luminosity to the stellar bolometric luminosity. We calculate lower limits on the amount 
of dust present in each system, after accounting for the effects of radiation pressure and 
Poynting-Robertson drag. Our analysis of the excess in the HD 17925 system points to the 
fact that it may be the star with the smallest amount of circumstellar dust detected among 
Vega-like systems to date. 

We present strong evidence for the existence of a previously undetected debris disk 
around HD 107146, and possibly one around HD 88638. Our detection of an excess in 
HD 107146 from archival IRAS FSC data suggests that there may be more stars with un- 
known debris disks in the deeper and higher signal-to-noise FSC, than what is currently 
known from analyses of the PSC. A methodical search for these in combination with 2MASS 
near-IR photometry, as already presented in Fajardo-Acosta, Beichman, & Cutri (2000) for 
the PSC, but watchful of the caveats addressed here and the ones described in Song et al. 
(2002), could therefore prove fruitful. Given that very few disks are found around late-type 
stars due to their intrinsic faintness, the results of such a search could significantly increase 
the number of known debris disks around solar- and later-type stars. Ground-based 10/iin 
and 20fim follow-up observations will be indispensable in discerning IRAS/ISO false posi- 
tives from bona-fide mid-IR excesses, as well as potentially successful in finding new excess 
candidates (e.g., as in the case of HD 88638). 

Based on our limited data, we discuss the possibility of breaking the degeneracy between 
debris disk temperature and mass in the fd vs. stellar age relation. We find that the power 
law index in the inferred dust mass M dust -age relation is steeper than in the /d-age case, and 
hence tentatively propose the use of M dust -age relation as a better tracer of the evolution 
of the amount of debris around main-sequence stars. However, a deficiency of the proposed 
power-law is that it does not take into account the absence of debris disk around stars which 
would be expected to harbor such based on their age. We argue that other factors, such as 
stellar multiplicity and spectral type, may also play a role. 
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Future space-based IR observations with SIRTF and SOFIA are expected to greatly 
expand the number of known Vega-like stars, by being sensitive to debris disks as faint as 
the zodiacal dust (fd ~ 8 x 10~ 8 Backman & Paresce 1993) around nearby stars as evolved 
as our Sun. These will provide the database necessary to answer the prominent questions 
about the processes and time-scales of evolution of debris disks around main-sequence stars, 
and will place our own solar system in the context of stellar and planetary evolution. 

We would like to thank Randy Campbell and Gregory Wirth for their assistance with 
LWS, Tom Hayward and Rick Burress for help with SC-10, and Jonathan Foster for help 
with the LWS data reduction. This research has made use of the IPAC InfraRed Science 
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in this paper were analyzed at ESA's ISO Data Centre at Vilspa, Spain. The publication 
makes use of data products from the Two Micron All Sky Survey, which is a joint project 
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Fig. la. — Spectral energy distributions of the observed stars. Our 10/zm data are plotted 
with solid circles, while data from the literature are plotted with open squares (with error- 
bars where available). Upper limits are plotted as downward facing triangles. Photospheric 
models (thin lines; photospheric temperatures listed under the object names) are from Ku- 
rucz (1979, 1992) and from Hauschildt et al. (1999). Dashed lines represent blackbody fits 
to the observed excess emission. The total flux densities (photosphere + blackbody) are 
plotted with thick lines. 
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Fig. 2.— Fits to the SED of HD 17925, where 250 K and 70 K blackbodies have been 
allowed to independently account for the excess. Symbols are as in Figure 1. The NextGen 
photosphere has been normalized to the J, H, and K band measurements, and the temper- 
ature and normalization of the blackbody fit to the excess have been allowed to vary freely. 
The two cumulative curves (thick continuous lines) fit the 10-60yum data comparatively 
well (Xpdf = 2-3 and 2.5, respectively). However, the hotter (250 K) blackbody produces a 
higher fractional excess fd than the cooler (70 K) one by half an order of magnitude, whereas 
the amount of dust Md ust required is less by nearly an order of magnitude. The best fit is 
achieved with a 140 K blackbody, corresponding to M dust 2 x 10~ 6 M e . 
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Fig. 3. — Top panel: f d vs. stellar age plot. Bottom panel: M disk vs. stellar age plot. In both 
panels solid dots represent stars with detected IR excesses, and are labeled. Down-pointing 
open triangles represent upper limits on fd, or on the dust mass Md us t, in systems with 
no excess IR emission, and are not labeled. Vertical errorbars are given for HD 17925 and 
HD 88638, corresponding to the range of inferred dust masses. The errorbars on the other 
points are ~ 0.5 dex in the top panel, and ~ 1.0 dex in the bottom panel. The data are 
from Tables 3 and 4, assuming greybody (GB) particle properties. 
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Table 1. Observational Epochs 
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2002 Jan 1 


SC-10 


field star 


RE J0137+18B 


RX J0137.6+1835B 


2002 Jan 1 


SC-10 


field star 


HD 60737 


HIP 37170 


2002 Jan 1 


SC-10 


field star 


HD 70573 


SAO 116694 


2002 Jan 3 


SC-10 


field star 


HD 70516 


HIP 41184 


2002 Jan 3 


SC-10 


field star 


HD 77407 


HIP 44458 


2002 Jan 3 


SC-10 


field star 


HD 88638 


HIP 50180 


2002 Jan 3 


SC-10 


field star 


HD 107146 


HIP 60074 


2002 Jan 3 


SC-10 


field star 



tCandidate TWA member (Sterzik et al. 1999). 



Table 2. Infrared photometry and properties of the observed sources 



Source 


J - H 


H -K a 


K s -N 


N 


AN 


*r s -[i2] 


sp. type 


age 


parallax 


Notes* 




(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 




(Myr) 


(mas) 




HD 145519 


0.18 ±0.04 


0.14 ±0.05 


-0.44 ±0.10 


7.04 ±0.10 


-0.10 ±0.13 


< 1.29 


B9V 


5-8 


6.3 ± 0.8 


1,2,3 


HD 109573A 


0.00 ±0.08 


0.00 ±0.08 


0.22 ±0.12 


5.58 ±0.10 


0.29 ±0.11 


< 1.25 


A0V 


10 ±5 


14.91 ±0.75 


1,4 


HD 147809 


0.22 ±0.06 


0.16 ±0.06 


-0.15 ± 0.12 


7.21 ±0.12 


-0.11 ±0.13 


1.91 ±0.11 


A1V 


5-8 


6.3 ± 0.8 


2,3 


HD 102647 


0.34 ±0.02 


0.01 ±0.01 


-0.08 ± 0.08 


2.06 ± 0.08 


-0.03 ±0.09 


0.03 ± 0.06 


A3V 


10-300 


90.16 ±0.89 


5,6 


CD -38°6968 


0.29 ±0.06 


0.03 ±0.05 


-0.12 ± 0.19 


8.60 ±0.19 


-0.10 ±0.19 


< 2.96 


G3 


< 100 


7.90 ±23.80 


7 


HD 107146 


0.27 ±0.03 


0.07 ±0.03 


0.19 ±0.09 


5.35 ±0.09 


0.11 ±0.09 


0.05 ±0.13 


G2V 


50-250 


35.07 ±0.88 


8,9 


HD 60737 


0.28 ±0.03 


0.06 ± 0.03 


-0.11 ± 0.11 


6.36 ±0.11 


-0.18 ±0.11 


< 0.73 


GO 


50-250 


26.13 ± 1.04 


9,10 


HD 70573 


0.28 ±0.04 


0.09 ± 0.04 


0.05 ±0.18 


7.14 ±0.18 


-0.01 ±0.18 


< 1.67 


G1/2V 


50-250 


11.30 ± 10.70 


9,11 


HD 70516 


28 + 03 


n 09 + 02 


KJ • i-'-t _1_ \J . VJ ,J 


fi 00 + 09 


(1 (17 + (1 (N 


n + 1 1 


GO 


1 00-300 


27 1 + 2 1 8 


9 


HD 77407 


0.27 ±0.03 


0.09 ±0.03 


-0.01 ± 0.08 


5.45 ±0.08 


-0.08 ±0.08 


0.07 ±0.11 


GO 


20-150 


33.24 ±0.91 


9,12 


HD 88638 


0.27 ±0.03 


0.11 ±0.03 


0.29 ±0.07 


6.04 ±0.07 


0.19 ±0.08 


< 0.81 


G5 


250-1000 


26.66 ±3.09 


9 


HD 17925 


0.39 ±0.02 


0.09 ±0.01 


0.02 ±0.10 


4.00 ±0.10 


0.14±0.11 


0.01 ±0.05 


K1V 


50-250 


96.33 ±0.77 


5,8,9 


HD 216803 


0.54 ±0.02 


0.11 ±0.01 


0.03 ±0.14 


3.80 ±0.14 


0.06 ±0.14 


0.07 ±0.06 


K4V 


200±100 


130.94 ± 0.92 


5,9,13 


CD -40° 8434 


0.73 ± 0.05 


0.28 ±0.05 


3.05 ±0.10 


5.49 ±0.10 


3.14 ±0.12 


3.56 ±0.09 


K5 


<io 




14 


RE J0137+18A 


0.62 ±0.03 


0.14 ±0.03 


0.26 ±0.14 


7.18 ±0.14 


0.15 ±0.15 


< 1.92 


K3V 


30-55 


15.8 ±5.0 


9,12,15,16,17 


RE J0137+18B 


0.61 ±0.03 


0.15 ±0.03 


0.32 ±0.14 


7.23 ±0.14 


0.22 ±0.15 


< 1.91 


K3V 


30-55 


15.8 ±5.0 


9,12,15,16,17 


AP 93 


0.50 ±0.03 


0.14 ±0.03 


< 1.18 


> 8.18 


< 1.1 


< 3.84 


K2 


80 


5.3 


18,19 


AA Tau 


0.87 ±0.05 


0.51 ±0.05 


3.23 ±0.08 


4.81 ±0.07 


3.07 ±0.09 


3.58 ±0.11 


MOV 


1 


7.1 


20 


DI Tau 


0.72 ±0.03 


0.21 ±0.03 


0.10 ±0.22 


8.29 ±0.22 


-0.08 ±0.22 




M0.5V 


0.6 


7.1 


21 


CD -33° 7795 


0.68 ±0.04 


0.24 ±0.04 


0.24 ±0.10 


6.50 ±0.10 


-0.11 ±0.10 


< 1.23 


M3V 


10 ±5 


18 ±3 


4,22 


Hen 3 -600 A 


0.62 


0.30 


3.04 ±0.20 


4.26 ±0.19 


2.54 ±0.10 


3.48 ± 0.09 


M4 


10 ±5 


18 ±3 


4,22,23,24 



tl. JHK photometry from Jura et al. (1993); 2. Assumed age for Upper Sco OB association is 5-8 Myr (de Geus et al. 1989); 3. Mean distance for Upper Sco 
from de Geus et al. (1989); Jones (1970); 4. Assumed TW Hya age is 10 ± 5 Myr; 5. JHK photometry from Aumann & Probst (1991); 6. Maximum age estimate 
from Song et al. (2001). Minimum age estimated from the lack of known accretion signature in this star; 7. Age estimate from Sterzik et al. (1999); 8. Age deduced 
from Li measurement of Wichmann, Schmitt, & Hubrig (2003); 9. Age estimate from comparison of unpublished Li A6707 equivalent width to open cluster sequence; 
10. Spectral type from Houck & Swift (1999); 11. Age deduced from Li measurement reported in Jeffries (1995); 12. Age estimate from Montes et al. (2001); 13. 
Age estimate from Barrado y Navascues (1998); 14. Age estimate based on Li measurement and spectral type estimate in Gregorio-Hetem & Hetem (2002); 15. The 
binary is unresolved in 2MASS. JHK S photometry from adaptive optics observations by Mctchcv & Hillenbrand, in prep.; 16. Spectral type from Jeffries (1995); 17. 
Distance estimate from Montes et al. (2001); consistent with being a nearby member of the Cas-Tau OB association (E. Mamajck 2003, private communication); 18. 
3cr upper limit at 10/im; 19. Spectral type deduced from effective temperature listed in Randich ct al. (1998); 20. Age estimate from Beckwith et al. (1990); 21. Age 
estimate from Meyer ct al. (1997); 22. Distance obtained as the average of the distances to the 4 TWA members with Hipparcos astrometry: TW Hya, HR 4796A, 
HD 98800, and TWA 9 (56, 67, 47, and 50 pc, respectively); 23. The binary is unresolved in 2MASS. JHK photometry from Webb ct al. (1999); 24. Spectral type 
from Torres et al. (2000). 
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Table 3. Fractional disk excesses 



Star Age [Myr] T dust [K] f d /d(Lit) Reference Excess? 



HD 145519 


5-8 


30 - 140 


2.7 - 8.0 x 10~ 4 




1 


Yes 


HD 109573A 


10 ±5 


110 ± 10 


6.1 x 10~ 3 


5 - 10 x 10" 3 


2,3 


Yes 


HD 147809 


5-8 


(300) t 


< 5 x 10~ 4 






No 


HD 102647 


10 - 300 


100 ± 10 


2.0 x 10~ 5 


1.2 - 1.9 x 10~ 5 


4,5 


Yes 


CD -38° 6968 


< 100 (300) 


< 10" 3 








No 


HD 107146 


50 - 250 


60 ± 10 


1.5 x 10- 3 






Yes 


HD 60737 


50 - 250 


(300) 


< 10~ 3 






No 


HD 70573 


50 - 250 


(300) 


< 10" 3 






No 


HD 70516 


100 - 300 


(300) 


< 10" 3 






No 


HD 77407 


20 - 150 


(300) 


< 10~ 3 






No 


HD 88638 


250 - 1000 


130 - 1500 


1 x 10~ 4 - 6 x 10~ 2 






Yes? 


HD 17925 


50 - 250 


30 - 250 


4 x 10~ 6 - 3 x 10~ 4 


1.3 x 10~ 4 


6 


Yes 


HD 216803 


200 ± 100 


(300) 


< 5 x 10~ 4 






No 


CD -40° 8434 


<io 


270 ± 20, 45 ± 5 


0.24 


0.34 


7 


Yes 


RE J0137+18A 


30 - 55 


(300) 


< 10" 3 






No 


RE J0137+18B 


30 - 55 


(300) 


< 10~ 3 






No 


AP 93 


80 


(300) 


< 7 x 10~ 3 






No 


AA Tau 


1 


1000 ± 100, 190 ±10, 40 ±5 


0.51 


0.31 


8,9 


Yes 


DI Tau 


0.6 


(300) 


< 3 x 10~ 3 






No 


CD -33° 7795 


10 ±5 


(300) 


< 10~ 3 






No 


Hen 3-600A 


10 ±5 


260 ±10, 30 ±10 


0.16 


0.21 


3 


Yes 



^For the stars without detected excesses, the listed upper limits of arc for assumed 300 K debris disks. 

References. — 1. Backman ct al. (1997), M. Meyer et al., in prep.; 2. Jura (1991); 3. Weinberger et al. (2002); 4. Jayawardhana et 
al. (2001); 5. Spanglcr ct al. (2001); 6. Habing et al. (2001); 7. Grcgorio-Hetcm & Hetem (2002); 8. Beckwith et al. (1990); 9. Chiang 
ct al. (2001). 



Table 4. Derived disk properties 



Star 


ft m j n 

H 


(a) 

H 




D [AU] 






Minimum M^ ust 


[M e ]t 




t PR [Myr] 




Minimum Mpe 

[M®]t 




BB 


GB 


ISM 


BB 


GB 


ISM 


BB 


GB 


ISM 


HD 145519 


5.7 


9.4 


700-30 


2500-50 


10 5 -10 3 


0.5-0.001 


6.5-0.003 


4000-0.08 


130-0.23 


1600-0.6 


10 6 -10 2 


0.009-0.027 


HD 109573A 


3.3 


5.5 


37 


93 


1100 


0.14 


0.83 


120 


0.38 


2.4 


340 


2.7 


HD 102647 


1.8 


3.1 


30 


100 


840 


1.6xl0~ 4 


1.7xl0~ 3 


0.12 


0.32 


3.5 


250 


0.025 


HD 107146 


0.36 


0.60 


28 


290 


890 


2.0xl0~ 3 


0.22 


2.0 


0.48 


51 


480 


0.85 


HD 88638 


0.22 


0.36 


4-0.03 


40-0.08 


70-0.09 


2xl0~ 6 -6xl0- 8 


2xl0~ 4 -4xl0 


~ 7 5xl0" 4 -5xl0" 7 


io- 2 -io~ 6 


l-io- 6 


3-10~ 5 


0.09-50 


HD 17925 


0.10 


0.16 


50-0.8 


1000-7 


2000-7 


5xl0- 6 -9xl0~ 8 


5xl0~ 3 -7xl0 


" 6 Ixl0" 2 -7xl0" 6 


2.0-0.0005 


2000-0.04 


4000-0.04 


2xl0~ 4 -2xl0- 2 


CD -40° 8434 






>0.85 


>6.9 


















AA Taut 






>0.077 


>9.1 






1.8xl0 5 












Hen 3-600A 






>0.96 


>6.3 



















tThe assumed mass of the Earth is M e = 5.9736 X 10 27 gram. 

tTotal dust+gas mass from Beckwith et al. (1990); Chiang et al. (2001). 
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Table 5. Assumed stellar parameters 



Star 


Adopted T cS 










(K) 


(Re) 


(Le) 


(M ) 


HD 145519 1 


11000 


2.2 


63 


2.6 


HD 109573A 


10000 


2.0 


34 


2.4 


HD 102647 


8500 


1.8 


15 


1.9 


HD 107146 


6000 


1.2 


1.7 


1.1 


HD 88638 


5500 


1.0 


0.88 


0.94 


HD 17925 


5000 


1.1 


0.37 


0.89 


CD -40°8434 2 


4400 


1.39 


0.64 




AA Tau 3 


4000 


2.1 


0.98 


0.67 


Hen 3-600A 4 


3200 


2.8 


0.72 





1 T e g and L* from de Geus et al. (1989). T c g is rounded 
to the nearest 500 K for compatibility with Kurucz 
(1979) models. 

2 Stellar parameters from Gregorio-Hetem & Hetem 
(2002). 

3 Stellar parameters from Beckwith et al. (1990). 

4 Stellar parameters estimated assuming bolometric 
correction of -3.24 (for M4 spectral type; Flower 1996), 
and F=12.04 (Torres et al. 2000). 



